Large bubble rupture sparks fast liquid jet 
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This Letter presents the novel experimental observation of long and narrow jets shooting out in 
disconnecting large elongated bubbles. We investigate this phenomenon by carrying out experiments 
with various viscosities, surface tensions, densities and nozzle radii. We propose a universal scaling 
law for the jet velocity, which unexpectedly involves the bubble height to the power 3/2. This 
anomalous exponent suggests an energy focusing phenomenon. We demonstrate experimentally 
that this focusing is purely gravity- driven and independent of the pinch-off singularity. 



Disrupting bubbles often exhibit violent jets during 
their ultimate fast out-of-equilibrium dynamics. Exam- 
ples include bursting sea-[Tll2] or champagne-bubbles [3j, 
collapsing cavitation clouds f4] or submarine explosions 
[5j up to the astronomically sized buoyant bubble in el- 
liptical galaxy M87 [6]. Among these violent jets, some 
are the signature of a finite-time singularity [71 occurring 
only for particular values of the control parameters [8j, 
but in many cases, including those described here, jet 
formation is simply the result of a nonsingular relaxation 
process 

Before bursting however, bubbles are not known to ex- 
hibit interface deformation as intense as the ones just 
depicted. Indeed, irrespective of their initial distortion 
p.0] small bubbles gently relax toward robust equilib- 
rium shapes in the form of spheres or oblate spheroids 
pT] . As regards the final form of larger bubbles, ini- 
tial shape matters. According to their initial appear- 
ance, weakly deformed bubbles smoothly evolve either 
to spherical caps [12j or toroidal bubbles fi3\ [T4], as the 
ones typically expelled by dolphins or divers [15]. Much 
less is known about initially strongly deformed bubbles, 
as their dynamics and outcome have barely been studied. 
Such bubbles are yet commonplace in bubbling systems 
typically found in geophysics or industry (e.g. occurring 
in glassy [16] or metallic [17J melts). The primary aim 
of the present study is therefore to investigate the tran- 
sient behavior of such large and highly deformed bubbles, 
which will prove to be as intense as bursting events. 

In this Letter, we report a surprising violent jet dy- 
namics following large bubble disconnection (Fig. [T]) in 
a bubbling experiment. The thin and concentrated jet 
developing inside the bubble possibly gives rise to liquid 
projections shooting out way above the free surface. We 
analyze experimentally this phenomenon and we provide 
evidences of the role played by gravity in the jet forma- 
tion. 

Our experiments consist in releasing gas bubbles from 
a submerged orifice in a viscous liquid. The liquid is 
contained in a transparent tank taken sufficiently large 
(20 cm X 20 cm x 25 cm) to rule out confinement ef- 
fects. It was observed that liquid slugs in the feeding line 
could make the bubble detachment frequency strongly 



fluctuate. Those slugs originate in a liquid invasion, ei- 
ther caused by capillarity or gravity. Capillary invasion 
is circumvented by the use of non-wetting aqueous liq- 
uids instead of oils. Therefore the liquids used in this 
study include sugar cane syrup (Canadou) of viscosity 
/i = 110 mPa.s and surface tension 7 = 90 mN.m~^ and 
three water-glycerol mixtures of viscosity fi = 140, 280, 
420 mPa.s and surface tension 7 = 65 mN.m~^. In order 
to both allow the development of the bubbles and mit- 
igate gravity-driven invasion (high hydrostatic pressure 
on the orifice), the height h of fluid above the injector is 
typically kept in the range 5-10 cm. Air is fed through a 
nozzle of diameter 1.8 mm. Injection is controlled by a 
mass flow meter (Alicat Scientific) that provides a con- 
stant flowrate Q by adjusting the air pressure. It allows 
us to achieve a wide range of airflows (from 0.01 to 10 
£.mm~^). The bubble and jet dynamics are analyzed 
through ultra- fast imagery. To do so, the tank is back lit 
and images are obtained at 4000 frames per second using 
a digital high-speed camera (Photron SA-5). 

Figure [T] illustrates a typical jetting event following the 
disconnection of a large gas bubble. The entire sequence 
lasts 18 ms and the bubble height, right after detach- 
ment, is H = 4.22 cm. Bubble pinch-off occurs on the 
first image and instantly a narrow, high-speed vertical 
jet shoots out inside the bubble. It is noteworthy that 
the bubble shape hardly deforms throughout the devel- 
opment of the jet, except near the fast recoiling conical 
rear. As the jet velocity Vjet = 2.69 m.s~^ is much higher 
than the bubble front velocity V^ront = 0-67 m.s~^, the 
jet tip reaches the top of the bubble in a few milliseconds. 
Here, the jet is suddenly stopped as it collides with the 
bubble wall, but in other cases the jet literally perforates 
the bubble and makes headway in the bulk liquid. 

In addition to the classical forces determining bubble 
volume [18j, the detachment of the bubble is here po- 
tentially influenced by several other complex processes, 
including viscous stresses associated with the large scale 
convection pattern, gas pressure at the inlet, wake of the 
foregoing bubble. . . However we argue in the following 
that detachment is primarily due to the gravity-driven 
collapse of the bubble neck. We derive a simple model 
for the bubble formation, considering principally the ex- 



FIG. 1. Time sequence of the jet developing inside a large bubble just after its release from a submerged nozzle. The large 
bubble pinches off at the first image and the back-to-equilibrium dynamics exhibits an intense and concentrated jet. The time 
lapse between the snapshots is At = 3 ms. The bubble height, right after detachment, is H = 4.22 cm, the bubble front velocity 
is Vfront = 0.67 m.s~^ and the jet velocity is Vjet = 2.69 m.s~^. The liquid viscosity is /x = 420 mPa.s. The gas is injected 
through an injector of diameter d = 1.8 mm at an airflow rate of Q = 4.4 ^.min~^. The fluid height is /i = 10 cm. 



panding and contracting liquid motions induced by the 
bubble growth and seal process. The ambient liquid is 
mimicked by a finite set of slices in pure radial motion 
around a pressurized cavity [19, 20j. Assuming a poten- 
tial flow evolution, we can obtain an equation for the 
position of the free boundary R{t) for each slice: 

V / Roo 2 p pR 

This is the 2D version of Rayleigh-Plesset equation de- 
scribing the liquid motion around a hollow cylinder 
[21] [22] . with p the liquid density, g the gravity, A = 
1 — R^/R^^ a confinement factor, Pbubbie = —P9H{t)-\-jhz 
the bubble pressure (taking the reference level at z = 0), 
Hi = 2/i?front the frout curvature, R^o a distant cut-off, 
and H{t) = l^ront t and z denoting respectively the height 
of the bubble and the local position of the slice. The 
bubble front is approximated by a frozen spherical cap 
of angle 20 rising at constant velocity. The cavity for- 
mation results from the air entrained behind it. The 
initialization procedure for each slice below this travel- 
ing bubble front simulates the deflection of fluid particles 
past it. At the trailing edge of the sphere portion, i.e. 
at the point R{t) = i^front sin 6 the following velocity is 
imposed: R{t) = dR/dt = Vcot^. Figure [i] shows a 
typical face-to-face comparison between the experiment 
and the model. Remarkably, the overall detachment se- 
quence is well captured even if the axial motions are dis- 



regarded [23]. This agreement between the experiment 
and a rough model containing as only driving force hy- 
drostatic pressure constitutes a strong evidence that the 
seal mechanism is gravity- driven. 

We now turn to the jet formation per se, occurring 
right after the cavity closing. Fig. [3^ shows the experi- 
mental relation between the jet velocity V^et and the bub- 
ble height right after pinch-off instant. Each circle 
(o) corresponds to a given experiment conducted with 
fixed ambient viscosity p = 420 mPa.s, injector diam- 
eter d = 1.8 mm and fluid height h = 10 cm (same 
physical parameters as in Fig. [T]). For all experiments, 
it was found from spatiotemporal diagram analysis that 
the velocity of the jet rapidly reaches a constant value 
after a transient regime. A clear trend appears from 
this figure: the jet velocity T/jet increases with the bub- 
ble height iJ, consistently with a gravity-powered mech- 
anism. Though exhibiting a visible tendency, the data 
are still quite spread out as for one value of the jet ve- 
locity correspond several bubble heights. To understand 
this spreading, let us examine the initial shape of some 
typical bubbles for a specific value Vjet — 2.7 m.s~^. In 
Fig.jsja) are represented three bubbles, marked from 1 to 
3, corresponding to the highlighted circles (bubble num- 
ber 3 is the same as in Fig. [T]). We observe that the 
three bubbles have a completely different shape, present- 
ing a wide range of front curvature. Interestingly, the 
largest bubble (3) also bears the highest front curvature. 
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FIG. 2. Comparison of experimental and simulated time 
sequences of bubble formation. In the 2D Rayleigh-Plesset 
model, the frozen spherical cap is materialized by a thick red 
line. The time lapse between each snapshot in the model and 
in the experiment is At = 11 ms. The physical parame- 
ters are the air flowrate Q — 3.80 ^.min~^, liquid viscosity 
ji — 280 mPa.s, density p — 1250 kg.m~^ and surface tension 
7 = 65 mN.m~^. The parameters of the model i^front = 4.18 
mm and Vfront = 0.785 m.s~^ are given by the experimental 
data. Roo = 18.3 mm and ^ = 1.1 rad are adjusted to the 
dynamics of the experiment. The slight shape dissimilarity is 
related to the absence of vertical motions in the theoretical 
model. Nonetheless, the seal timescale and overall dynamics 
are in fair agreement. 



while the smallest bubble (1) is the less curved. Yet in 
spite of these differences their jet velocities are the same. 
This geometrical competition between bubble height and 
front radius of curvature suggests a physical competi- 
tion between an hydrostatic driving force and a capillary 
quenching effect: the main source of pressure difference 
between the liquid and the bubble is gravity, and capil- 
larity tends to lessen this difference by pressurizing the 
bubble. In order to take both effects into account, we 
represent in Fig. 3]^a) the jet velocity V-^^t versus the cap- 
illary corrected bubble height = H — 2^/ pgR^ronu 

is a natural lengthscale in the sense that the bubble 
pressure is approximately —pgH^ (see orange circles ( ) 
related to the bottom axis). Upon using instead of 

the three tracked experiments now get together as 
indicated by the arrows in Fig. |3ja). A similar gather- 
ing is observed for all experiments. This indicates the 
relevance of the parameter and points to a pressure- 
driven mechanism for the jet formation, as is common in 
impact-driven [9j or gravity-driven jets [24J. 
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FIG. 3. (a) Jet tip velocity (Vjet) as a function of heights: 
H (o) height of the bubble after detachment (top axis) and 

— H — 2j/ pgRhont ( ) taking into account the capillary 
pressure term (bottom axis). The data spread as a function 
of H significantly decreases when plotted as a function of 

. The three bubbles presented in the bottom right cor- 
ner labelled 1, 2 and 3 have the same jet velocity, but for a 
different height H (highlighted in blue). When plotted with 
respect to , they all gather. The viscosity of the liquid is 
p = 420 mPa.s. (b) Vjet as a function of for four different 
viscosities /i : 110 mPa.s ( ), 140 mPa.s ( ), 280 mPa.s (4), 
420 mPa.s ( ). The surface tension and density are 7 = 65 
mN.m~^ and p = 1250 kg.m~^ except for the lowest viscosity 
( ) for which 7 = 90 mN.m"^ and p = 1350 kg.m"^. The 
injector diameter is d = 1.8 mm and the fluid height is /i = 10 
cm. 



The influence of viscosity p and surface tension 7 in the 
development of the jet has been investigated as well by 
conducting similar experiments in various liquids. The 
results are summarized in Fig. [sjb). For each viscosity 
the same rise in jet velocity V^et with is observed, al- 
beit with a general shift to higher with more viscous 
liquids. It is worth noting that jets still exist for values 
of the control parameters corresponding to the shaded 
area in Fig. [sj^b). But in this region, the jet velocity 
never reaches a constant value. Such unsteady jets were 
systematically disregarded during postprocessing. This 
shaded area corresponds to a range of 'capillary Froude 
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number' V-^et/ below a typical value of 5. There, 
the feeding velocity ]^et is too close to the retraction ve- 
locity ~ Y^7/pr, with r the jet radius, hence the un- 
steadiness. 

Now taking the natural gravito-inertial velocity \fgH^ 
and the nozzle radius Rq as relevant scales of the prob- 
lem, we non-dimensionalize our results. The rescaled jet 
velocity or Froude number Fr = V-^et/y/gH^ is plotted 
against the dimensionless height / Rq in the bottom 
right inset of Fig.[4]for all viscosities, surface tensions and 
nozzle radii. Surprisingly enough, the six corresponding 
curves increase with respect to the dimensionless height. 
Actually we might have expected from a balance between 
buoyancy and inertia to observe a constant Froude num- 
ber [24J. Rather we observe a linear variation of the 
Froude number with /Rq^ identical for all sets of ex- 
periments. Moreover a shift to lower non-dimensional 
velocities is noticeable as the viscosity is increased. This 
points to a functional dependence of the Froude number 
with / Ro and the dimensionless viscosity as follows: 



jet 



Ro 
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pVgH^Ro 



with a a non-dimensional constant. Analysis of the ex- 
perimental data reveals that the form of function {x) 
is actually a linear function with /3 constant. This 
allows us to rewrite the preceding relation with a single 
identical offset for all viscosities: 

Re = a Ar — /3, 

where we have introduced the jet Reynolds number Re = 
pVjetRo/ p and Archimedes number Ar = py/gH^H^ / ji. 

Figure [4] represents the relation between the jet 
Reynolds number and Archimedes number for all the con- 
ducted experiments. The collapse of the six series of ex- 
periments corresponding to five different viscosities, two 
surface tensions and densities and three different nozzle 
diameters is excellent. The relation between the jet ve- 
locity and the capillary corrected height of the bubble for 
the whole range of physical parameters is finally captured 
with the simple law Re = a Ar - /3, taking a = 0.33 and 
/3 = 13. In systems exhibiting a balance between buoy- 
ancy and inertia, the potential volume energy ~ pgH 
is converted into kinetic volume energy ~ pV'^ and the 
velocity typically scales as H^^'^: large Taylor bubbles 
in tubes jEj, gravity waves in shallow fluid layers [25j, 
inertial gravity current [26j . . . Surprisingly, the here ob- 
tained scaling law shows a dependence of the jet velocity 
V-^et with H^^^^ rather than H^^^^ . This anomalous ex- 
ponent suggests an energy focusing phenomenon whose 
source is now discussed. 

Does the pinch-off singularity play a role in this focus- 
ing? In order to investigate the role of detachment, we 
have designed and carried out a new model experiment, 
free of pinch-off. This experiment consists in blowing air 
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FIG. 4. Jet Reynolds number Re = pVjetRo/l^ versus 
Archimedes number Ar = p^JgH* / ii for different viscosi- 
ties, surface tensions, densities and nozzle radii. The signifi- 
cance of the symbols is the same as in Fig. [sj^b) and two series 
were added corresponding to two different injector diameters : 
2.15 mm [<) and 1.19 mm ( ) with viscosity, surface tension 
and density respectively fi — 320 mPa.s, 7 = 65 mN.m~^ and 
p — 1250 kg.m~^. The equation of the oblique dashed line 
is Re = a Ar - /3, with a — 0.33 and (3 — 13. The bottom 
right inset shows the Froude number Fr = V-^et / \fgH^ versus 
the nondimensional bubble height /Rq. The lines act as 
eye guides and present the same slope a. The top left inset 
presents the results of a different experiment where the jet 
is formed due to the gravitational collapse of a free surface 
hollow. By taking H the depth of the cavity and L the maxi- 
mum cavity diameter as relevant scales, jet Reynolds number 
Re = pVjetL/p is plotted versus Archimedes number Ar = 
py/gHH/p for two different viscosities p : 220 mPa.s ( ) and 
426 mPa.s ( ). The two snapshots display a typical cavity 
(H — 5.61 cm and L — 12.6 cm) and the resulting jet. 



over a free liquid surface so as to form a depression, mim- 
icking the conical rear of the bubble. Upon the collapse of 
this hollow, a gravity- driven jet develops (see snapshots 
fig.[4|. But there again, the same scaling for the jet veloc- 
ity with iJ^/^ is observed (top left inset). Note also that, 
in sharp contrast with singular focusing behaviors [7J, 
the cavity reversal observed here is not unlike the shell 
eversion process in that jet eruption/curvature reversal 
happens before collapse. These rule out both pinch-off 
[27] and curvature singularities [7J in the observed energy 
focusing and demonstrate that it is merely a signature of 
the gravitational cavity collapse. 

In conclusion, the present experimental work reports 
on the violent dynamics exhibited during the relaxation 
of an initially large oblate bubble. The intense and nar- 
row jets developing inside the detaching bubbles follow a 
surprising dependence with H^^^^ . The related focusing 
of energy is not a consequence of the detachment singu- 
larity, as proven by the persistence of the scaling in a 
pinch-off-free setting. Consequently, in this experiment 
at least, the jet does not keep the footprint of the singu- 
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larity as in e.g. Zeff et al. [7], Bartolo et al. f5S] or Gekle 
and Gordillo [29j. Instead, global conservation rules take 
over the jet dynamics. The scaling law also suggests the 
existence range of this type of jet, with a threshold in 
Archimedes given by j3/a {c^ 40 in our case). Below 
this value, the deformation of the bubble rather tends to 
the unsteady liquid tongue of Walters and Davidson [13j. 
Alternatively, the threshold defines a critical viscosity- 
dependent bubble height necessary to observe those liq- 
uid jets. Finally, in this paper we have focused on the 
jet development before its collision with the bubble front. 
An open remaining question is whether the jet could have 
a sufficient energy to perforate the bubble and reach the 
free surface. In our experiments we observed that strong 
liquid projections can appear above the free surface for 
large airflow rate. We are currently investigating the 
highly unclear link between the jet presented here and 
these liquid projections. 
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